HAMILTONIAN FORMALISM FOR SOLVING THE VLASOV-POISSON 
EQUATIONS AND ITS APPLICATION TO THE COHERENT 
BEAM-BEAM INTERACTION 
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Abstract 



A Hamiltonian approach to the solution of the Vlasov- 
Poisson equations has been developed. Based on a nonlin- 
ear canonical transformation, the rapidly oscillating terms 
in the original Hamiltonian are transformed away, yielding 
a new Hamiltonian that contains slowly varying terms only. 
The formalism has been applied to the coherent beam-beam 
interaction, and a stationary solution to the transformed 
Vlasov equation has been obtained. 



1 INTRODUCTION 

The evolution of charged particle beams in accelerators 
and storage rings can often be described by the Vlasov- 
Maxwell equations. At high energies the discrete-particle 
collision term [Jj]] comprises a small correction to the dy- 
namics and can be neglected. Radiation effects at suffi- 
ciently high energies for leptons can be a significant fea- 
ture of the dynamics, and should be included in the model 
under consideration. 

The Vlasov-Maxwell equations constitute a considerable 
simplification in the description of charged particle beam 
propagation. Nonetheless there are only a few cases that are 
tractable analytically. Therefore, it is of utmost the impor- 
tance to develop a systematic perturbation approach, able 
to provide satisfactory results in a wide variety of cases of 
physical interest. 

Particle beams are subject to external forces that are of- 
ten rapidly oscillating, such as quadrupolar focusing forces, 
RF fields, etc. In addition, the collective self-field excita- 
tions can be rapidly oscillating as well. A typical example 
is a colliding-beam storage ring device, where the evolu- 
tion of each beam is strongly affected by the electromag- 
netic force produced by the counter-propagating beam. The 
beam-beam kick each beam experiences is localized only in 
a small region around the interaction point, and is periodic 
with a period of one turn. 

In this and other important cases one is primarily inter- 
ested in the long-time behavior of the beam, thus discard- 
ing the fast processes on time scales of order the period 
of the rapid oscillations. To extract the relevant informa- 
tion, an efficient method of averaging is developed in the 
next section. Unlike the standard canonical perturbation 
technique |Q, the approach used here is carried out in 
a "mixed" phase space (old coordinates and new canonical 
momenta), which is simpler and more efficient in a compu- 
tational sense. The canonical perturbation method devel- 



oped here is further applied to the coherent beam-beam in- 
teraction, and a coupled set of nonlinear integral equations 
for the equilibrium beam densities has been derived. 

2 THE HAMILTONIAN FORMALISM 

We consider a A^-dimensional dynamical system, de- 
scribed by the canonical conjugate pair of vector variables 
(q, p) with components 



q 
p 
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The Vlasov equation for the distribution function /(q, p; t) 
can be expressed as 



df 
Of 



where 



_ dF_ dG_ _ dF dG 

dqi dpi dpi dqi 



(2.2) 



(2.3) 



is the Poisson bracket, H (q, p; t) is the Hamiltonian of the 
system, and summation over repeated indices is implied. 
Next we define a canonical transformation via the generat- 
ing function of the second type according to 



S = S(q,P;t), 
and assume that the Jacobian matrix 

d 2 s 



J*,(q,P;£) = 



dqidPj 



is non-degenerate with 



det (Jirf^O, 



(2.4) 



(2.5) 



(2.6) 



so that the inverse J i4 1 exists. Let us also define the distri- 
bution function in terms of the new coordinates (Q, P) and 
the mixed pair (q, P) as 



/(q,p;t) = /o(Q,P;«) = F (q,P;t). 



(2.7) 



The new canonical variables (Q, P) are defined by the 
canonical transformation as 



Pi 



dS_ 

dqi ' 



OS 



(2.8) 



Because 

dpi 



8 2 S 



dPj dqidPj 



dp 

dpj 



zii = 7- 



(2.9) 



we can express the Poisson bracket in terms of the mixed 
variables in the form 



where the new Hamiltonian K is a function of the new 
canonical pair (Q, P), such that 



1 ' J J 3* \dqidPj dP 3 dq t 



(2.10) }C(W P S,P;t) = ^+H(ci,V q S;t), 



Differentiation of Eq. (2.8) with respect to time t, keeping 
the old variables (q, p) fixed, yields 



d 2 S d 2 S fdPj 



dqidt dqidPj \ dt 



0, 



dQi 
dt 



d 2 S fdP, 



d 2 S 



dPidt dpdP 3 V dt , 



9P± 
dt 



= -J< 



_! d 2 S 



j i 



dq^t 



(2.11) 
(2.12) 
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Our goal is to express the Vlasov equation (2.2) in terms of 
the mixed variables (q, P). Taking into account the identi- 
ties 



dQ t 
d qj 



d 2 S 
dq-jdPi 



— Jji 



wr J ^ (2 - l4) 



dfr ^dFo 
dQi Jij By' 

dfo dF df d 2 S 



dp t dp dQj dP.dPj ' 



(2.15) 
(2.16) 



we obtain 



Of 



dt 



dfo , dfofdQ 



dt dQ, V dt 



dfo (dp 
dPi \ dt 



<n> 



dFp xf dFp d 2 S dFp d 2 S 
dt Jl \ dq t dtdP 3 dPj dtdqt 

dS~ 



= dFo 
dt 

Furthermore, using the relation 



F , 



dt 



(2.17) 



[f,H] = [F ,ni 



(2.18) 
(2.19) 



where 

W(q,P;t) = fl-(q,V,5;t), 

we express the Vlasov equation in terms of the mixed vari 
ables according to 



where 



dF Q 
dt 



%P;i) 



[^0,^=0, 



(2.20) 



dS_ 
~dt 



H(q,y q S;t) (2.21) 



is the new Hamiltonian. For the distribution function 
/o(Q, P; t) depending on the new canonical variables, we 
clearly obtain 

dfo 



^ + [/o,/C] = 0, 



(2.22) 



(2.23) 



and the Poisson bracket entering Eq. (2.22) has the same 



form as Eq. (2.3 ), expressed in the new canonical variables. 



3 COHERENT BEAM-BEAM 
INTERACTION 

As an application of the formalism developed in the pre- 
vious section we study here the evolution of two counter- 
propagating beams, nonlinearly coupled by the electromag- 
netic interaction between the beams at collision. For sim- 
plicity, we consider one-dimensional motion in the vertical 
(q) direction, described by the nonlinear Vlasov-Poisson 
equations 

dfk 
dO 



[fk,H k ]=0, 



(3.1) 



dq 2 



47r dpf 3 -k(q,p; 0), 



(3.2) 



where 



Hk = y (p 2 + q 2 ) + \ k 6 p (6)V k (q; 9) (3.3) 



is the Hamiltonian. The notation in Eqs. ( |3.l[ ) - (3.3) is 
the same as in Ref. [Bh. Our goal is to determine a canon- 
ical transformation such that the new Hamiltonian is time- 
independent. As a c onseq uence, the stationary solution of 
the Vlasov equation ( Z.20| ) is expressed as a function of the 
new Hamiltonian. Following the procedure outlined in the 
preceding section we transform Eqs. (3.1) - (3.3) according 
to 

0, (3.4) 



Fp {k \lC k 



dS k 

de 

d 2 V k 



eH k [q,^;6)=JC k (q,P), (3.5) 



4tt dP 



<}!Sk Ft k \q,P), (3.6) 



dqdP~° 

where e is formally a small parameter, which will be set 
equal to unity at the end of the calculation. The next step is 
to expand the quantities S k , K, k and V k in a power series in 
e as 



S k = qP + eG k 1] 

IC k = e/Cj^ - 

V k = V k + eV k {1) - 

where _ 

d 2 V k 
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(3.7) 

(3.8) 
(3.9) 

(3.10) 



Subs titution of the above expansions (3.7) - (3.9) into Eqs. 
(3.5) and (3.6) yields perturbation equations that can be 
solved successively order by order. The results are: 



First Order. 0(e) 



4"(9.P) = y(P 



-•■« + £&(.), 



V k W (q;9)^0. 



Second Order: 0(e 2 ) 



Kf> (q,P)=Q, 



G k 2) (q,P; 



2tt 



n#0 



Afc^ T 7( 2 ) 



And 



where 



g> 2 T4 (2) 
dq 2 



*KV£{q)ldPFt k \<l,P) 



(3.11) 
(3.12) 
(3.13) 

(3.14) 
(3.15) 

(3.16) 
(3.17) 



Third Order. 0(e 3 ) In third order we are interested in 
the new Hamiltonian, which is of the form 



4 3) (<?^) 



C(2) ^' 2 (g)-214 (2) (g) 



where ((z) is Riemann's zeta-function 



(3.18) 



(3.19) 



n=l 



4 THE EQUILIBRIUM DISTRIBUTION 
FUNCTION 

Since the new Hamiltonian Kk is time-independent (by 

(k) 

construct ion) , the equilibrium distribution function Fq 
[see Eq. (3.4 )] is a function of the new Hamiltonian 



F^ k \q,P)=g k (K, k ), 



(4.1) 



where 



K k (q,P) = ^(P*+ q 2 ) 



At 
2tt 



V k (q) 



Xl *<{2)\vf{q)-2Vl%j 



Air 2 



(4.2) 



Integrating Eq. 1 [ > over P we obtain a nonlinear integral 
equation of Haissinski type [|[] for the equilibrium beam 



(k) 

density profile 



g ( k) (q)= ldPg k {K k ), 



(4.3) 



where 



IC k (q,P) = ^(P 2 +q 2 



q-q 



\at k) tf) 



2A^C(2)^ fc (g) > (4.4) 



^(<z) = y dq'dq"Z{q-q',q' - q") 



,(3-A), 



x^o V)<??~*V)> (4-5) 
Z(u,v) = sgn(u)sgn(v) — 2|u|<5(w). (4.6) 

5 CONCLUDING REMARKS 

We have developed a systematic canonical perturbation 
approach that removes rapidly oscillating terms in Hamil- 
tonians of quite general form. The essential feature of this 
approach is the use of mixed canonical variables. For this 
purpose the Vlasov-Poisson equations are transformed to 
mixed canonical variables, and an appropriate perturbation 
scheme is chosen to obtain the equilibrium phase space 
density. It is worthwhile to note that the perturbation ex- 
pansion outlined in the preceding section can be carried 
out to arbitrary order, although higher-order calculations 
become very tedious. 

The canonical perturbation technique has been applied to 
study the one-dimensional beam-beam interaction. In par- 
ticular, rapidly oscillating terms due to the periodic beam- 
beam kicks have been averaged away, yielding a time- 
independent new Hamiltonian. Furthermore, the equilib- 
rium distribution functions have been obtained as a general 
function of the new Hamiltonian, and a coupled set of inte- 
gral equations for the beam densities has been derived. 
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